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Abstract 

Based on a Boltzmann-like traffic equation for aggressive drivers we construct in this paper a second- 



< ^ ■ order continuum traffic model which is similar to the Navier-Stokes equations for viscous fluids by applying 
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two well-known methods of gas-kinetic theory, namely: the Chapman-Enskog method and the method of 
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O moments of Grad. The viscosity coefficient appearing in our macroscopic traffic model is not introduced in 
an ad hoc way - as in other second-order traffic flow models - but comes into play through the derivation of a 

£h ■ first-order constitutive relation for the traffic pressure. Numerical simulations show that our Navier-Stokes- 
like traffic model satisfies the anisotropy condition and produces numerical results which are consistent with 
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our daily experiences in real traffic. 
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I. INTRODUCTION 



Because of the great variety of phenomena present in the motion of vehicles along highways or in 
urban roads, traffic dynamics has attracted the attention of a large number of researchers in the 
last decades. Traditionally, three different approaches can be used to study traffic flow problems, 
namely: a purely microscopic approach in which the acceleration of a driver- vehicle unit is determined 
by other vehicles moving in traffic flow, a macroscopic approach which describes the collective motion 
of the vehicles as the one- dimensional compressible flow of a fluid, and a mesoscopic approach which 
specifies the individual behavior of the vehicles by means of probability distribution functions. 

A reading of relevant literature shows that since 1955, when Lighthill and Whitham [l| proposed 
the first continuum model to describe traffic flow, much progress has been made in the development 
of macroscopic (fluid-type) models on the one hand, and of microsocopic (follow-the-leader) models 
on the other hand. The first mesoscopic (or gas-kinetic) traffic flow model just appeared in 1960, 
when Prigogine and Andrews Q wrote a Boltzmann-like equation to describe the time evolution of a 
one-vehicle distribution function in a phase-space where the position and the velocity of the vehicles 
plays a role. Until the 1990s, mesoscopic traffic models do not get much attention from scientists 
due to their lack of ability to describe traffic operations outside of the free-flow regime. Addition- 
ally, compared to macroscopic traffic flow models, gas-kinetic traffic models have a large number of 
independent variables that increase the computational complexity. However, in the last decade, the 
scientific community's interest by mesoscopic traffic models resurrected with the publication of some 
works that apply these models to derive macroscopic traffic models (see, for example, the papers 
of Helbing 3|, Hoogendoorn and Bovy 4], Wagner et al. Q]). Macroscopic equations for relevant 
traffic variables can be derived from a Boltzamnn-like traffic equation by averaging over the instan- 
taneous velocity of the vehicles. This is a well-known procedure in the kinetic theory, nevertheless 
its application leads to a closure problem, i.e., there are some quantities which must be evaluated 
with constitutive relations in order to obtain a system of closed equations. The analogy with well- 
established methods of the kinetic theory of gases - such as the Chapman-Enskog method jf| or 
the method of moments of Grad [3] - gives us a clue to proceed, provided we have at least a local 
equilibrium distribution function. 

By assuming that motorists drive aggressively 8| we apply both the Chapman-Enskog method and 
the method of moments of Grad to construct a second-order continuum traffic model which is similar 
to the Navier-Stokes model for viscous fluids. In the Chapman-Enskog method, constitutive relations 



are constructed at sucessive levels of approximation by expanding the distribution function in power 
of the mean free path, while in Grad's moment method the distribution function is approximated 
by an expansion in orthonormal polynomials where the coefficients are related to the moments of 
the distribution function. The dependence of our constitutive relation for the traffic pressure with 
the velocity gradient in non-equilibrium situations, drive us to define a traffic viscosity coefficient 
which, in our case, depends on the traffic state through the density and the mean velocity of the 
vehicles. As several second-order continuum traffic models, there exists in our Navier-Stokes-like 
traffic model a characteristic speed that is greater than the average flow velocity. The existence 
of this faster characteristic speed means that the motion of the vehicles will be influenced by the 
traffic conditions behind them. This seems to be a drawback of our second-order model since one 
fundamental principle of traffic flow is that vehicles are anisotropic and respond only to frontal stimuli. 



101 ] . this faster characteristic 



However, as pointed out by Helbing and Johansson j^] and J. Yi et al. 
speed does not represent a theoretical inconsistency to our Navier-Stokes traffic model since it is 
related to an eigenmode that decays quickly and, therefore, it cannot emerge by itself. Besides, 
we check the anisotropic behavior of our second-order continuum traffic model by performing the 
simulations of a traffic situation where a discontinuity is present, namely: the removal of a blockade 
scenario. 

The organization of the paper is as follows: in Section 2 we present our kinetic traffic model 
for aggressive drivers. In Section 3 we construct a second-order continuum traffic model - which is 
similar to the Navier-Stokes equations for viscous fluids - by applying both the Chapman- Enskog 
method and the method of moments of Grad. In Section 4 we present the results of our numerical 
simulation. Finally, we give in Section 5 a summary. 



II. KINETIC TRAFFIC MODEL 

In a mesoscopic description the state of a vehicle at a given instant of time is represented by points 
in phase space, where the coordinates are the position x and the velocity c of the vehicle. Thus, in 
analogy with the kinetic theory of gases, an one- vehicle distribution function f(x, c, t) can be defined 
in such a way that f(x, c, t) dx dc gives at time t the number of vehicles in the road interval between 
x and x + dx and in the velocity interval between c and c + dc. For an uni-directional single-lane 
road without entrances and exits, the one-vehicle distribution function satisfies the kinetic traffic 
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equation [11] 



dt dx dc V dt 
where the interaction term 
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Q(f,/) = J (l-p)(c'-c)/(i,c,«)/(i,c',«)dc'- / (l-p)(e-e')f(x,c,t)f(x,c',t)<k' (2) 

describes the decelaration processes due to slower vehicles which can cannot be immediately over- 
taken. The first part of the interaction term corresponds to situations where a vehicle with velocity 
d must decelerate to velocity c causing an increase of the one-vehicle distribution function, while 
the second one describes the decrease of the one-vehicle distribution function due to situations in 
which vehicles with velocity c must decelerate to even slower velocity d . The derivation of the in- 
teraction term is based on the following hypotheses: (i) vehicles are regarded as point-like objects, 
(ii) the slowing down process has the probability (1 — p), where p is the probability of passing, (iii) 
the velocity of the slow vehicle is not affected by interactions or by being passed, (iv) there is no 
braking time, (v) only two-vehicle interactions are considered and (vi) vehicular chaos is assumed, 
in such way that the two-vehicle distribution function can be factorized. The individual acceleration 
term appearing on the left-hand side of the kinetic traffic equation can be modeled by assuming that 
vehicles moving with velocity c accelerate exponentially to their desired velocity c = c (x,c,t) with 
a relaxation time r, i.e. 

dc c c . . 

It ~ ~- [6) 

In fact, the desired velocity of the vehicles is determined by the average balance among several traffic 
parameters like legal traffic regulations, weather conditions, road conditions and drivers personality, 
i.e., it is a phenomenological function. Despite of the variety of traffic parameters that determines 
the desired velocity of the vehicles, we shall consider in this work the simple relation 8] 

c = wc, (4) 

where w = w(x, t) > 1 is a positive parameter that depends on traffic conditions along the road. On 
driver's level, relation (4) indicates that the desired velocity of the vehicles increases as their velocity 
increases, which is a common feature of aggressive drivers. By assuming that the aggressiveness 
parameter w depends on position and time only through the vehicular density, we shall consider the 
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FIG. 1. Dependence of the aggressiveness parameter on the vehicular density 
following constitutive relation 

Po ~ Pc 

w = l + (w c -l)^(^P-) P< (5) 

Pc \P -PJ 

where p is the jam (or maximum) density, p c is the critical vehicle density (i.e. the density value 
which maximizes the aggressiveness parameter) and w c = w(p c ) is the maximal value for drivers' 
aggressiveness. Expression (5) tell us that: (i) when traffic is very dilute the aggressiveness parameter 
must tend to one, since in this limit vehicles move along the highway at its desired speed, (ii) as the 
vehicle density increases, the velocity of the vehicles starts to decrease, and drivers try to compensate 
this situation by increasing their aggressiveness, (iii) when the density exceeds its critical value, it 
becomes increasingly difficult for the vehicles to move at its desired speed and, in this case, the 
drivers have no choice except to reduce its aggressiveness, (iv) when the highway is fully occupied, 
the vehicles can no longer move and the aggressiveness parameter must tend to one. Figure 1 
shows the dependence of the aggressiveness parameter on the vehicular density for p c /p = 0.2 and 
w c = 1.2. Lastly, it is important to mention that the constitutive relation (5) is not a unique model 
to the aggressiveness parameter, but just a sound one. 
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III. SECOND-ORDER CONTINUUM MODEL 



The kinetic traffic equation (1) allows the derivation of balance equations for macroscopic traffic 
variables like the vehicular density 



o 



p(x,t) = I f(x,c,t)dc (6) 

and the average velocity 



oc 



c^f^dc. (7) 
p(x,t) 







The integration of the kinetic traffic equation over all values of the actual velocity of the vehicles 
yields the continuity equation 

dp dpv 



dt dx 

while the traffic momentum equation 



0, (8) 



dpv d , o N w — 1 

~k + m ( p + ^ = p ~ ~ v ~ p{ ~ p)w (9) 

follows through the multiplication of the kinetic traffic equation with c and the integration over all 
values of the actual velocity of the vehicles. In the derivation of the traffic momentum equation we 
have introduced the traffic pressure 

POO 

w(x,t)= / (c- v) 2 f(x,c,t)dc, (10) 
Jo 

and used relations (3) and (4). At this point it is important to emphasize that the balance equa- 
tions (8) and (9) can only be obtained if the one-vehicle distribution function satisfies the following 
boundary conditions: 

]imf(x, c, t) — and lim f(x, c, t) — 0. (11) 

c— >0 c— >oo 

Based on the continuity and momentum equations, we can construct a second-order continuum traffic 
model by specifying the traffic pressure in terms of the vehicular density, the average velocity and 
their spatial gradients. Since there are a variety of possible constitutive relations which can be 
borrowed from fluid dynamics, we shall restrict ourselves here to the derivation of a constitutive 
relation for the traffic pressure which is similar to the usual Navier-Stokes relation for ordinary 
viscous fluids, i.e., a constitutive relation written in terms of the density, the average velocity and 

6 



their first-order spatial gradients. One can achieve this goal by applying, for example, the Chapman- 
Enskog method or the method of moments of Grad, as they are developed in the kinetic theory of 
gases. In the Chapman-Enskog method, constitutive relations are constructed at sucessive levels of 
approximation by expanding the distribution function in power of the mean free path, while in Grad's 
moment method the gas-kinetic equation is replaced by a set of balance equations for the moments 
of the distribution function. To close this set of equations, the distribution function is approximated 
by an expansion in orthonormal polynomials, where the coefficients are related to the moments of 
the distribution function. Then, by applying an iteration procedure in the resulting system of field 
equations, it is possible to derive first-order constitutive relations. Here, we shall apply both methods 
to derive a first-order constitutive relation for the traffic pressure which is similar to the Navier-Stokes 
relation for viscous fluids. 

A. Chapman-Enskog Method 

The basic idea of the Chapman-Enskog method is to expand the distribution function into a power 
series in the rarefaction parameter e as 



where f {r) represents successive approximations to the distribution function. In traffic flow problems, 
the rarefaction parameter is a dimensionless number defined as 

<=\ (13) 

where L is the road length and A denotes the mean distance traveled by a vehicle between successive 
interactions. Based on the value of the rarefaction paramater, it is possible to identify two traffic 
flow regimes: (i) the rarefied regime, when £ < 1, and (ii) the continuum regime, when e ^> 1. In 
the rarefied regime, the processes of acceleration and decelaration of the vehicles can be neglected 
and the general solution of the kinetic traffic equation in this regime reads 




oo 



(12) 



r=0 



f(x,c,t) = f(x - ct,c,t = 0). 



(14) 



Otherwise, in the continuum regime, the processes of acceleration and deceleration of the vehicles can 
not be neglected, so that the distribution function is determined by both the individual acceleration 
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law and the interaction term. To determine the distribution function in this limit, it is appropriate 
to write the kinetic traffic equation (1) in the following dimensionless form 

% + C % + £ ¥c{ f Jt) =£ SI {l ' ?){C '~ C) °' ^ (15) 
where the position is given in units of the road length, the velocity is given in units of the average 
velocity, the time is given units of the mean flow time (the mean time required for vehicles to traverse 
the road) and the vehicular acceleration is given in units of v 2 /X. 

Integral equations for the distribution function approximations can be easily obtained if we insert 
the expansion (12) into the dimensionless kinetic traffic equation (15) and equate equal power of the 
scale parameter. Hence, we obtain 

§- c (r^) = QU (o \r) (i6) 

and 



+ + £[/ w ^l =Q(/ W ,/ W ) (r>l) (17) 



dt dx dc V dt 



where 



r POO 

Q(f°\ f (r) ) = V / (1 - P) (c - c) /(»>/'(-«> dc'. (18) 

By assuming that the vehicular density and the mean velocity are only determined by the zeroth- 
order approximation of the distribution function, the solution of the integral equation (16) with the 
individual acceleration law (3) leads to 

a £ /«p 
1 [a) v \ v / V v / 

where 

p (1 — p) vr , . 

a = — . 20 

w — 1 

Expression (19) for the zeroth-order distribution function tell us that the velocity of the vehicles 
is gamma-distributed with a shape parameter a and a rate parameter (3 = v/a when drivers drive 
agressively on the highway. To gain an insight into the shape parameter, let us calculate the velocity 
variance (or velocity dispersion) in the zeroth-order approximation, i.e., 

(C-V) 2 f {0) dc 2 

= -• (21) 

f m dc a 
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From the above expression we verify that the velocity variance depends quadratically on the mean 
velocity, a fact which can be used to identify the inverse of the shape parameter as the prefactor of the 
velocity variance. The experimental data reported by Shvetsov and Helbing 12] demonstrate that 
the prefactor of the velocity variance is almost constant at the low-density region, otherwise it can be 
taken as a function of the vehicular density In this work, we shall take this prefactor or, rather, the 
shape parameter as a constant whose value is greater than unity, so that our gas-kinetic-like traffic 
model is restricted to low densities. 

In order to determine the first-order approximation to the one-vehicle distribution function we 
introduce expression (19) into the integral equation 



df 



Of 



df m d 
dx dc 



f 



id 



dc 
dt 



Q(F\f w ) 



(22) 



which follows from (17) by taking r = 1. By utilizing the constraints 

poo poo 

/ f^dc = and / c f r) dc = 
Jo Jo 

which are valid for r > 1 we obtain 



(23) 



/ 



(0) 



p p\v 



) dx 



a 



1 



a 



1 



2 dv 
dx 



+ 



d_ 



f 



in 



dc 
dt 



-p(l-p)( c - w )/« (24) 



where the dot denotes the material time derivative. Integration of equation (24) over all values of 
the actual velocity leads to the continuity equation 



P 
P 



Or 
dx 



(25) 



Besides, if we multiply equation (24) with c and integrate over all values of the actual velocity we 
get 



v 
a — 

v 



v dp dv 
p dx dx 



a(l 



(26) 



If equations (25) and (26) are used to eliminate the material time derivatives appearing in equa- 
tion (24) we obtain 



/ 



(0) 



a[ - -1 

v 



2[--l]-l 



- + -(/ 

dx dc 



HI) 



w — 1 



T 



a[ - - 1 



w — 1 



T 



f 



(i) 



Fa(l-p)v(^-l) £°(^-l)V C1 W. 



(27) 
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Based on equation (27) we shall write the first-order approximation to one-vehicle distribution func- 
tion as 



= / 



(0) 



a + Oi - lj + a 2 - 1 j 



dv 
dx 



(28) 



where a , a 1 and a 2 are position and time-dependent coefficients to be determined. Insertion of 
expression (28) into constraints (23) and equation (27) leads to 



a 2 
a 



T 



(29) 



2 a 2(w-l)' 

Note that the dependence of the expansion coefficients with position and time comes through the 
parameter w that measures the average agressiveness of drivers. Hence, the first-order approximation 
to the distribution function reads 



/ 



in 



/ 



(0) 



at'- -I 2 



21 



dv 
dx 



(30) 



2 (w - 1) |_ \v 

Once the one-vehicle distribution function is a known function of the basic traffic variables and 
their first-order spatial gradients, the traffic pressure can be directly calculated as 



(c-vYU {0) + f W )dc 



,2 ™,2 
2— T 

a a 



pir n pir /a+l\dv 



a J dx 



(31) 



where 



(32) 



2(w-iy 

The constitutive relation (31) for the traffic pressure has a similar form to the Navier-Stokes relation 
for viscous fluids since in non-equilibrium situations both constitutive relations depend on the velocity 
gradient. Based on this similarity, it is possible to define a traffic viscosity coefficient 

,2 _ Ll . 



pv fa 

PL = fi[p, v) = 2 r 

a \ a 



(33) 



which depends on the traffic state through the vehicular density and the average velocity. At this 
point, it is important to mention that a similar constitutive relation for the traffic pressure was derived 
by Velasco and Marques Jr. 8| by applying a simplified version of the Chapman-Enskog method to 
the reduced Paveri- Font ana traffic equation [l^]. In their formalism, the collective relaxation time 
t appears as a free adjustable parameter of order of the mean vehicular interaction time, and it was 
introduced by means of a relaxation approximation performed in the interaction term. 
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B. Grad's Moment Method 



In Grad's moment method a macroscopic description of traffic flow is based on macroscopic traffic 
variables like the vehicular density, the average velocity and the central moments of the distribution 
function 



m k (x,t)= {c-v) k f(x,c,t)dc (k>2). (34) 
Jo 

The balance equations governing the dynamical behavior of these macroscopic traffic variables are 
the continuity equation (8), the traffic momentum equation (9) and the balance equations 

dm k d dv m k -i dw w - 1 

-7T- + tt {rn k v + m k+ x) + km k - k k m k 

at ox ox p ox t 

= -p(l - p) ( m k+1 - k^^w ). (35) 



P 

In deriving the balance equations (35) for the central moments we have multiplied the kinetic traffic 
equation (1) by (c — v) k , and subsequently integrated with respect to the actual velocity of the 
vehicles. 

Clearly, we can see that the balance equations (8), (9) and (35) form a non-closed system of field 
equations for the determination of the moments p, v and m k , since the balance equation for the 
central moment m k contains the central moment m k+ i which is not a priori related to the lower 
order moments. The dependence of the central moment m k+ i upon the moments p, v and m k can 
be attained if we know the distribution function as a function of p, v and m k . In the method 
of the moments such normal solution is found by an expansion around a zeroth-order distribution 
function (19), i.e., we write the distribution function as 

oo 

f(x,c,t) = r(x lC ,t)J2Cn(x,t)P n (c) (36) 

n=0 

where C n (x,t) are position and time-dependent expansion coefficients and P n (c) are orthonormal 
polynomials in the actual velocity of the vehicles. Since in the zeroth-order approximation the 
velocity of the vehicles is gamma-distributed, it is possible to construct the orthonormal polynomials 
Pn(c) by applying the condition 

poo 

<P(s)P n (s)P m (s)ds = 5 nm (37) 
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where s = ac/v is the dimensionless instantaneous velocity and 

„a—l p —s 

*(*) = -f^tt- (38) 
r(a) 

is the probability density function of the gamma distribution. From the orthonormality condition 
(37) we verify that the first polynomials read 

P (s) = 1, (39) 

Pi(s) = S ~^, (40) 



1 } y/2a(a + l) V ' 

n , , s 3 - 3 (a + 2)s 2 + 3 (a + l)(a + 2)s - a (a + l)(a + 2) , . 

P3 S = /r 7 \w ■ 42 

v/6a(a + l)(a + 2) 



We can easily verify from the above expressions that the orthonorma 
to the associated Laguerre polynomials (see the textbook of Arfken 
formula 



polynomials P n (s) are related 



141 ] ) and they are given by the 



P (s) = -tWLJ F ^ *L (s^e-s) (43 ) 
n[S) s «-i e -»V n\T(a + n)ds n[ j ' 1 j 

By using the orthonormality condition (37) the position and time-dependent expansion coefficients C n 
can be determined as follows: 

poo °° poo 

/ P n (c)f(x, c,t)dc = P J2 C m / $(s)P n (s)P m (s) ds = P C n . (44) 

m=0 Jo 

Thus, we conclude that the expansion coefficients C n are related directly to the moments of the 
distribution function. For example, the first coefficients read 

C = 1, (45) 

C t = 0, (46) 



fi=j / +1 ;- gt , (47) 

2(a + 1) w 



_ 2a f <p-<Po o w-w \ 

3 " A/ 3(a + l)(a + 2)^ O 6 w 
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where 



oc 



4>{x,t) = m 3 (x,t) = / (c-vyf(x,c,t)dc (49) 
Jo 

is the third-order central moment. Besides, 

zu (x,t)= [ (c-v) 2 f 0) (x,c,t)dc= P — (50) 
Jo at 

and 

O (M)= nc-vfr(x,c,t)dc = 2^ (51) 
Jo « 

are the values of the second and third-order central moments in the zeroth-order approximation. In- 
sertion of the position and time-dependent coefficients into the expansion of the distribution function 
allows us to write it in terms of the velocity polynomials and the macroscopic traffic variables. Note 
that each coefficient in the expansion of the distribution function introduces a new macroscopic traffic 
variable, so that it is possible to choose which relevant variables we want to use in our macroscopic 
traffic description. 

Let us now construct a continuum traffic flow model based only on three traffic variables, namely: 
the vehicular density, the average velocity and the traffic pressure. In this case, the balance equations 
governing the dynamical behavior of traffic variables are the continuity equation (8), the momentum 
equation (9) and traffic pressure equation 

dvo d , , . „ dv n w — 1 ,„ . , , 

~at + di {wv + <p)+ w Tx ~ 2 —r w = _p(1 " p)0 - (52) 

The balance equations (8), (9) and (52) becomes a system of field equations for the determination of 
p, v and w, if a relationship can be established between these variables and the third-order central 
moment 0. In order to achieve this goal, expansion (36) for the distribution function is taken with 
C n (x, t) = for n > 3, so that we have 

f _ f( o) L s 2 -2(a+ l)s + «(« + !) w-w \ 

1 1 I 2(a + l) w J ' 1 ) 

A comparison of the distribution function for the first order Chapmann-Enskog expansion (30) and 
the one obtained by Grad's method (53) shows that the last expression (i) involves an extended set of 
variables, note that in this formalism the traffic pressure is an independent variable, and (ii) depends 
only on the macroscopic variables not on the gradients. If the traffic pressure is calculated assuming 
(53), we realize that the result is different from (31). To reconcile these differences we will carry 
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out an standard iterative procedure known in the kinetic theory of gases as the Maxwellian iteration 
procedure. 

Insertion of the distribution function (53) into expression (49) leads, after a simple integration, 
to the following constitutive relation for the third-order central moment: 

= 3— (w- -w j . (54) 



w \ 3 

If we introduce the constitutive relation (54) into the balance equations (8), (9) and (52) we obtain 
a system of field equations for p, v and w or, equivalently, for p, v and w, where w = w — w is the 
so-called traffic pressure deviator. Hence, after some algebra, we get 

— + v—+ — = (55) 
dt dx dx ' 

dv w dp . w dv dw w — 1 . , . ^. . , 

" * + f £ + (a + 2) + & = » - " (1 - p)( ^ + ^ (56) 

dw fa + A\cj) Q dw rt fa + 2\^dv n fa + l\dv w , r _, 

3 )w— + 2w { J — — — . (57) 



dt \ 2 / w dx V a J dx V a J dx r„ 

As in the kinetic theory of gases, one can transform the balance equation (57) into an approximate 
constitutive relation for the traffic pressure deviator by applying a method akin with the Maxwellian 



iteration procedure 15j. For the first-iteration step we insert, on the left-hand side, the value of the 
traffic pressure deviator in the zeroth-order approximation, namely w = 0, and get, on the right-hand 
side, the first iterated value 

^ pv 2 (a. + 1\ dv , 

™=- 2 -<— )*■ < 58) 

The above expression for the traffic pressure deviator is identical to that one obtained by the 
Chapman-Enskog procedure, a fact that allow us to affirm that Chapman- Enskog expansion and 
Grad's moment method are both physically and mathematically equivalent, at least in first-order 
approximation. 

C. Navier-Stokes-like Traffic Model 

Insertion of the constitutive relation (31) for the traffic pressure into the balance equations (8) 
and (9) leads to the so-called Navier-Stokes-like traffic model which can be written in the following 
conservative form 
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where 



U 



( 



F(U) 



/ 



and S(U) 



\pv 2 + pc 2 y 



bv. 



(60) 



In the above expressions we have introduced the traffic sound speed c(v) = ^dvj Q /dp , the optimal 
velocity function u(p,v) = wv — r(l — p)ca () and the anticipation coefficient 

a — 1 \ <9zz7 



< 0. 



(61) 



2 J dv 

In contrast to others macroscopic traffic models, we see that our optimal velocity function does not 
depend only on the vehicular density, but also on the average velocity and that such dependence is 
explicitly determined by the average desired velocity of the vehicles reduced by a term arising from 
deceleration processes due to vehicle interactions. Besides, in our macroscopic traffic flow model, 
traffic viscosity is not introduced in an ad hoc way - as in some of the most popular second-order 
contiunuum traffic models - but it comes into play via an interation procedure and reflects the way 
drivers anticipate traffic situation on the basis of second-order spatial changes in the average velocity. 

Finally, we close this section by remarking that the eigenvalues A of the conservative flux-Jacobian 
matrix 

/ 1 \ 

(62) 



A(U) 



<9F(U) 



V 



a + 1 „ a 
v 1 2 — 
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determine how traffic information (such as slown-downs and speed-ups) is transmitted in a traffic 
stream. These eigenvalues, also known as characteristic speeds, are found by setting 

det|A(U) - AI| =0, (63) 

where I is the identity matrix. Hence, our Navier-Stokes-like traffic model has two real and distinct 
characteristic speeds, namely: 

A, 



a 



1 =f + l 



a 



(64) 



Note that one of the characteristic speeds is larger than the average traffic flow velocity indicating 
that traffic disturbances propagate in the downstream direction. This fact has been criticized by 
Daganzo as a major deficiency of most second-order traffic models. However, as shown by Helbing 
and Johanson [q| in a recent publication, this fact does not constitute a theoretical inconsistency of 
our second-order macroscopic traffic model since the perturbation that travels faster than traffic 
decays quickly. 
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IV. NUMERICAL EXAMPLE 



Since the Navier-Stokes-like traffic equations (59) form a hyperbolic system of partial differential 
equations which involve smooth as well as discontinuous solutions, it is suitable to use a conservative 
numerical scheme to solve them accurately. However, as pointed out in the literature, conservative 
numerical schemes are only appropriate for solving the homogeneous part of the system, i.e. when the 
source term S(U) vanishes. In order to deal with the source term, the following fractional-step 
scheme [17J 

At 
Ax 

At 



ur = u° 



(Fi+1/2 — Fi_i/ 2 ) 



(65) 



U 



n+l 



U- + — (S(Uf 



S(U*)) 



(66) 



can be employed, where Ax is the grid size, At is the time increment and Fi+1/2 = F(Uf, Uf +1 ) is the 
intercell numerical flux. The intercell numerical flux F i+1 / 2 approximates the time-integral average 
of flux across the cell interface between cells i and i + 1, where the i-th cell is given by the spatial 
interval between £i_i/ 2 and Xj+i/2- One can find in the literature several classes of methods that 
provide stable and consistent approximations to the intercell numerical flux. Here, we determine 
the numerical flux by employing a first order upwind-type scheme based on Roe's flux-difference 



splitting [18 



The application of Roe's scheme to the homogeneous version of our macroscopic traffic model 
yields a conservative method whose numerical flux function is computed as 

F(UQ + F(U i+1 ) 1^ 



Fi+1/2 



(67) 



where A k , <j k and e k are respectively the eigenvalues, the wave strengths and the right-eigenvectors 
of the flux-Jacobian matrix A(U) at Roe's average state U. For the Navier-Stokes-like traffic model 
these quantities are expressed as: 

a + 1 =F y/a + 1 



Ai, 



0\ 



=F 



a 



2 + 1 



a 

A(pv) a 
v 



■v. 



1 ± VoTI 



a 



Ap 



a 



1 =f Va+T _ 
v 



a 



(68) 
(69) 

(70) 
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FIG. 2. Time evolution of the vehicle density for the removal of blockade scenario 

where AU = U; + i — Ui. As usual, Roe's average state is determined by imposing the shock- capturing 
property 

AF = F(U i+1 ) - F(Ui) = A(U) AU. (71) 

From (71) we verify that the first equation is identically satisfied, while the second one gives the 
following average for the velocity: 



y/Pi + VPi+i 

Note that no representation for average traffic density is required in our Navier-Stokes-like traffic 
model, since Roe's average (72) for the flow velocity completely defines the flux-Jacobian matrix. 
Concerning now the discretization of the source term we take 

blUij-pi - - — + ^ , [76) 



where 



H{p u vi) + /i(pi±i,fi±i) ( . 

A*i±i/2 = = • (74) 
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Completing the description of the numerical scheme it is worth to mention that the numerical stability 
of explicit difference schemes is ensured by the Courant-Friedrich-Lewy condition 

^-max |A k | < 1 (75) 

AX k 

which guarantees that information propagates only through a single cell at each time step. 

We illustrate now the numerical scheme presented above by considering a traffic situation where 
a queue of nearly motionless vehicles is present in a certain road region. At the initial time, the 
blockade at the head of the queue is removed and vehicles flow into the empty part of the roadway. 
For simulation of this traffic scenario, we consider the following initial conditions on a 20 km circular 
road: 

V if 2.5km <,< 7.6km 

p(x, 0) = ^ and v(x,0) = v e (p(x,0)), (76) 

, elsewhere 

where p* is the vehicle density in the queue and v e (p) is the equilibrium speed-density function. Here, 



we assume an equilibrium speed-density relation of the form (see Bando et al. 191 ] ) 



v e (p) = (tanh (p /p - a) + tanh (a)) , (77) 

where v is the free-flow velocity and a is a positive dimensionless constant. Furthermore, we impose 
periodic boundary conditions and choose the following values for the model parameters: 

a = 125, f o = 30m/s, p = 0.2veh/m, p* = 0.198 veh/m, a = 3.9 and r = 8s. (78) 

Regarding the probability of passing it is usual to assume that this quantity depends only on the 
vehicular density in a linear way. However, as pointed out by Hoogendoorn and Bovy [2^] in their 
report, an expression for the probability of passing that depends both on the vehicular density and 
the mean velocity can be obtained if we set u(p, v) = v e (p). By equating the optimal velocity to the 
equilibrium speed-density function we are in fact replacing the individual (microscopic) process of 
deceleration by some collective (macroscopic) relaxation process to an equilibrium traffic state. 

Figures 2 and 3 show that after the removal of the blockade the vehicles at the head of the 
queue move downstream into the empty region with the free-flow velocity, while vehicles at the 
tail of the queue remain at their location. Although the traffic conditions upstream are free-flow 
we observe from these figures that vehicles do not flow backwards into the empty region, a fact 
that allow us to say that our Navier-Stokes-like traffic model satisfies the anisotropy condition and 
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FIG. 3. Time evolution of the traffic flow for the removal of blockade scenario 
produces numerical results which are similar to traffic operations in real-life traffic. Similar results 



were obtained by Hoogendoorn 21[ through numerical simulations performed on the single user-class 
version of his macroscopic multiple user-class traffic flow model. Hoogendoorn's macroscopic multiple 
user-class traffic model was derived from mesoscopic principles which encompass contributions of 
drivers acceleration towards their user-class specific desired velocity and contributions resulting from 
interactions between vehicles of the same and different classes. Besides, the velocity variance is 
introduced as an additional basic field describing deviations from the average velocity within the 
user- classes. 



V. SUMMARY 

By appling both Chapman-Enskog expansion and Grad's moment method we constructed a second- 
order continuum traffic model which is very similar to the Navier-Stokes model for viscous fluids. The 
constitutive relation for the traffic pressure obtained by the method of moments of Grad is identical to 
that one obtained by the Chapman-Enskog procedure, a fact that allow us to affirm that Chapman- 
Enskog expansion and Grad's moment method are both physically and mathematically equivalent, 
at least in first-order approximation. Besides, in contrast to others second-order macroscoscopic 
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traffic models, our traffic viscosity coefficient - which depends on the traffic state through the vehicle 
density and the mean velocity - is not introduced in an ad hoc way, but comes into play through the 
derivation of a constitutive relation to the traffic pressure. Numerical simulation of a traffic scenario 
where a discontinuity is present show that our macroscopic traffic model satisfies the anisotropy 
condition and produces numerical results which are similar to traffic operations in real-life traffic. 
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